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Abstract 



The incident energy at which the azimuthal distributions in semi-central heavy 
ion collisions change from in-plane to out-of-plane enhancement - E tra n is studied as 
a function of mass of emitted particles, their transverse momentum and centrality 
for Au+Au collisions. The analysis is performed in a reference frame rotated with 
the sidewards flow angle (@fi ow ) relative to the beam axis. 

A systematic decrease of E tran as function of mass of the reaction products, their 
transverse momentum and collision centrality is evidenced. 

The predictions of a microscopic transport model (IQMD) are compared with the 
experimental results. 

Key words: NUCLEAR REACTIONS, E=90-400A MeV; semi-central collisions; 
flow angle, azimuthal distributions, transition energy; Quantum Molecular 
Dynamics model; nuclear matter Equation of State 



PACS: 25.75.Ld;25.70.Pq 



1 Introduction 

A decade ago, two types of azimuthal anisotropics have been evidenced in 
heavy ion collisions [1-4]. Following their shape relative to the reaction plane 
and the main mechanism behind them they have been baptized as in-plane 
and out-of-plane enhancement or rotational-like [1] and squeeze-out [2-4] phe- 
nomena, respectively. How the transition from one to the other scenario takes 
place as a function of incident energy, for different collision geometries and 
types of the reaction products was a question addressed few years later for 
a light [5] and heavy system [6]. Soon, it was realized that a detailed study 
of such a transition from in-plane to out-of-plane emission as a function of 
incident energy could give more insight into the relative contribution of the 
attractive and repulsive forces, the lifetime of the emitting source, its rota- 
tional energy and expansion dynamics. Therefore many studies concentrated 
on this subject [7-14]. 

Below this transition energy, E tran , the experimental data show an in-plane en- 
hancement of the azimuthal distribution. Theoretical calculations established 
that a rotating compound system created by the mean field, dissipating angu- 

1 Corresponding author: GSI, Planckstr. 1, D-64291 Darmstadt, Germany; 
Email: A.Andronic@gsi.de; Phone: +49 615971 2769; Fax: +49 615971 2989 
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lar momenta and excitation energy via particle emission, stays at the origin 
of this observation [7,15]. 

Increasing the beam energy the mean field which contributes to the formation 
of a rotating compound system becomes less important and processes based 
on nucleon-nucleon interaction start to be predominant. The hydrodynamic 
approaches used to predict different phenomena which take place in the rela- 
tivistic regime of heavy ion collisions have shown that a preferential emission 
of the "squeezed" participant zone in the free phase space, not hindered by the 
spectator nuclear matter, takes place. The out-of-plane emission of the nuclear 
matter has been called squeeze-out [16] and has been confirmed experimentally 
by Plastic Ball [2,4] and Diogene [3] Collaborations few years later. Afterwards, 
the generality of this phenomenon was evidenced experimentally. Thus, pions 
[17,18], kaons [19], neutrons [20,21], light particles [22,23], proton-like parti- 
cles [8] and intermediate mass fragments [11,12] show a similar squeeze-out 
pattern in the azimuthal distributions. Detailed theoretical investigations of 
this phenomenon were also undertaken, using microscopic transport models 
[24-30]. 

The existence of an azimuthally symmetric flow evidenced in central heavy 
ion collisions [32-35] can be regarded as the extreme case of the out-of-plane 
flow. In other words, the squeeze-out is the result of the expansion of the hot 
and compressed participant zone in the presence of the projectile and target 
spectators. Although it is obvious that the cleanest signal on the collective 
expansion of hot and compressed baryonic matter can be obtained from central 
collisions, one has to accept that this is an extreme case and the study of 
a rotating hot and compressed object remains an appealing subject. If the 
angular momenta and the shadowing realized by the spectator matter which 
come into the game bring more complications at first glance, they can be used 
as internal clocks to gain more information on the expansion dynamics. 

The incident energy where the effects of these three competing processes, 
expansion, rotation and shadowing, compensate each other could be looked as 
a benchmark observable for deeper understanding of the expansion dynamics. 

In this paper we present the results of a multidimensional analysis of the 
transition energy as a function of the mass of emitted particles, their transverse 
momentum and for different collision centralities. Early FOPI studies [6,11,12], 
done using the experimental data from Phase I of the detector, have been 
continued with the data obtained using the full phase space coverage of Phase 
II of this device [9,37,38]. Relative to the Phase I FOPI studies, the present 
data make possible the azimuthal distributions studies in a reference frame 
with the z axis along the sidewards flow direction, with no upper limitation on 
transverse momentum and as function of the mass of light charged particles. 
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Chapter 2 describes in more details the experiment and the experimental data 
analysis. The experimental results are presented in chapter 3. Chapter 4 is 
dedicated to the comparison of the results with the predictions of a microscopic 
transport model. Conclusions are presented in chapter 5. 



2 Experimental Details 

2.1 Setup 

The present experimental configuration [39,40] was used to study Au+Au, 
Xe+Csl and Ni+Ni collisions at 90, 120, 150, 250 and 400-A MeV, aiming 
to continue our studies started with Phase I of the FOPI [36] detector with 
much better phase space coverage and as a function of the baryonic number 
involved in the process. The beam energies are the mid-target energies, taking 
into account the upstream energy loss of the beam. The results presented in 
this paper refer to the Au+Au system. 

Since details of the first phase can be found in our earlier publications, es- 
pecially in ref. [36], we describe here the main features of the Central Drift 
Chamber (CDC), the main extension of Phase I of the FOPI facility [39,40] 
used in this experiment. The CDC is a drift chamber that performs the track- 
ing of the path followed by all charged reaction products emitted in the polar 
angular range 33° < 9 lab <150°. The path being curved in the magnetic field 
of a superconducting solenoid, one can use the transverse momentum and spe- 
cific energy loss for mass identification of the reaction products. Full azimuthal 
coverage is realized by this tracking device. The detector is subdivided in 16 
sectors, each containing 60 sense wires and 60 potential wires, all aligned par- 
allel to the beam axis. The sense wires are resistive and readout on both ends, 
which gives the possibility to reconstruct the position along the wire of the 
hit via charge division. The reconstruction of the initial position of the ionized 
cell along the ionizing particle path is obtained making use of the electrons 
drift velocity of about 43.7 /im/ns and of the Lorentz angle a L =32°. As the 
direction of the drift electrons cannot be measured, each hit has a symmet- 
ric partner relative to the sense wires plane (mirror hit). They combine and 
form mirror tracks. In order to have a criterion to distinguish between real 
and mirror tracks, the sense wires planes are tilted 8° relative to the radial 
direction starting from the beam axis. Within this geometry the mirror tracks 
do not originate anymore from the target. Furthermore, the sense wires within 
a plane are alternatively staggered by ±100 /im, which could be also used to 
reject the mirrored tracks. The individual signal wires of 50 /im thick NiCr- 
based alloy with a resistance of 500 Q/m vary between 86 and 190 cm. The 
mixture used for the chamber consists of 88% Ar, 10% isobuthane and 2% 
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CH 4 at a slight overpressure. 



The CDC was operated in a "split-mode" in order to increase its dynamical 
range. While the most outer 30 potential wires were operated at a nominal 
voltage of -1.55 kV, optimum for a good resolution for energetic light particles 
(Z=l,2), the most inner 30 potential wires had a lower voltage of -1.1 kV, such 
that the amplitude of the signals corresponding to highly ionizing fragments 
(Z=3-6) to be within the dynamical range of our sampling convertors, Flash 
ADCs (8 bit, non-linear). The relative momentum resolution u Pt /p t varies from 
4% for p t < 0.5 GeV/c to about 12% for pt=1 GeV/c. A helium bag placed 
between the target and the forward subdetector of Phase I [36], following the 
conical shape of the CDC forward endcap, was used in order to decrease the 
energy thresholds for the intermediate mass fragments flying in the forward 
direction at 6^ <30°. It is worth to mention here that within the phase space 
covered by the Phase I subdetector the reaction products are identified by 
their charge using specific energy loss and time-of-flight information, while the 
fragments detected by the CDC are identified by their mass, as the information 
on the specific energy loss and magnetic rigidity was available. 



Au target 




0.2 | 0.4 

90 1 50 250 400 



0.6 



0.2 | 0.4 

90 1 50 250 400 



Fig. 1. FOPI Phase II acceptance as a function of the transverse momentum per 
nucleon (pt/A) and laboratory reference frame rapidity for Z=l and Z=3 reaction 
products. The cm. rapidities corresponding to different incident energies are marked 
by arrows. 

The phase space coverage of the whole device used in the present experiment, 
in a transverse momentum - p t versus rapidity - y representation, can be 
followed in Fig. 1. The borders between different subdetectors at 1.2°, 7°, 30° 
and 33° are represented by thick lines. The shadowed area between 30° and 33° 
corresponds to the region between the Plastic Wall (PW) and CDC which is 
not covered at all. The low thresholds (LT) are represented by thin continuous 
lines. They correspond to the measured values (dotted lines) after correction 
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for the energy loss in the media before the active detectors. This was done 
assuming that the interaction has taken place in the middle of the target. For 
polar angles smaller than 30° the low threshold corresponds to the momenta 
for which the particles reach the PW in order to deliver a time signal while 
the specific energy loss is obtained from a layer of ionization chambers or thin 
plastic scintillators in front of the PW [36]. For polar angles larger than 30° LT 
was considered to correspond to the momenta for which the ionizing particles 
could reach the outer radius of the CDC. 

The dashed lines correspond, for the forward detector, to the momenta at 
which the plastic scintilators are just penetrated by the corresponding par- 
ticles. In this case the energy loss correction is negligible. For the CDC the 
dashed lines correspond to the momenta at which the energy loss corrections 
are below 1%. 



2.2 Analysis techniques 

The combination of the very high hit densities into tracks associated to the 
true particles is a serious challenge of a tracking algorithm [40,41]. The data 
analyzed in the present paper have been obtained using for the CDC a local 
tracking method - "track-following" - in a version developed within the FOPI 
collaboration [41]. 

Fig. 2 shows a representation where the performance of the algorithm in track 
resolution can be followed. It corresponds to a projection in a plane perpen- 
dicular to the beam axis of the trajectories associated to the charged particles 
seen by the CDC in a typical event from central Au+Au collisions at 250- 
AMeV. The straight lines of dots are representing the anode (sense) wires, 
while the rest of symbols are the registered hits. Grey symbols are belonging 
to reconstructed tracks and the open squares are their mirrored counterparts. 
The dark full squares are the leftovers remained unassigned to any track. 

The mass identification performance of CDC using specific energy loss 
(dE/dx) versus transverse momentum divided by charge (pt/q) of the detected 
particles can be followed in Fig. 3. A cross-check of the present calibration of 
CDC [42] was done using comparison with our earlier data measured with Si- 
Csl telescopes [43]. Good agreement between the two sets of data, in spectral 
shape and yields, was obtained. However, using only the CDC information, 
the A=3 branch contains a mixed contribution from 3 He and tritium (t) frag- 
ments (as can be seen in Fig. 3). Consequently, the momentum associated to 
3 He fragments will be underestimated by a factor of two if q—1, corresponding 
to hydrogen isotopes, is considered. Based on the phase space distribution of 
3 He and t obtained in one of our previous experiments [43], the contribution 
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Fig. 2. CDC xy plane cross section for a Au+Au, E=250-A MeV central event. 

of 3 He to A=3 yield for a given pf^ range (using q—1) is smaller that 15%. As 
it can be seen later, this amount of contamination does not influence different 
systematics of A=3 fragments. They nicely fall between those corresponding 
to A=2 and A=4 fragments. 

Most of the events were registered under "medium bias" (MB) trigger, corre- 
sponding to different values for the outer Plastic Wall multiplicity as a function 
of energy and mass of the colliding systems, such to reduce the background 
contribution written on tape, estimated from target free runs to be below 5%. 

Fig. 4 shows the multiplicity distribution of the reaction products detected 
and identified by CDC, CMUL, as function of Erat = Y^i^±,%/ J2i^\\,i (the 
sums run over all products detected in an event) for the energy of 250-A MeV 
under hardware "medium bias" trigger. As expected, besides a clear corre- 
lation between these two observables, producing the ridge of the displayed 
distribution, strong fluctuations are visible. The trends observed in this rep- 
resentation are in good qualitative agreement with the microscopic transport 
model prediction, particularly the IQMD version [44], extensively used for the 
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Fig. 4. CMVL-Erat correlation for 
Au+Au, E=250-AMeV. 



Fig. 3. The mass identification perfor- 
mance of CDC. The scale on z is loga- 
rithmic. 

present analysis. At large impact parameters (low CMUL and Erat values) the 
particle multiplicity has a higher selectivity of the collision geometry, while to- 
wards higher centrality (larger CMUL and Erat values) the Erat does better. 
Following similar recipe as the one used in Phase I for PMUL [45], a CMUL 
selection was devised in order to select the collision geometry. 

Au+Au E=250 AMeV 



-Q 
—1 




-Q 



60 80 
CMUL 




Erat 



Fig. 5. CMUL and Erat distributions for 250-A MeV. Different centrality bins used 
in the present analysis can be followed. 

In Fig. 5 we show the CMUL distribution (both for minimum and medium 
bias triggers) and the Erat distribution for CM4+CM5 bins. The CM2 and 
CM3 regions have been used in order to select impact parameters in the range 
of 6-8 fm and 4-6 fm, respectively and ER4 and ER5 for 2-4 fm and 0-2 
fm, respectively. The impact parameter values have been obtained using a 
geometrical "sharp cut-off" approximation for the reaction cross section. 
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Au+Au E=150 AMeV IQMD HM 




2 4 6 8 10 2 4 6 8 10 2 4 6 8 10 



Impact Parameter - b (fm) 

Fig. 6. The impact parameter distributions extracted from the IQMD (HM) model 
at the energy of 150-A MeV for the experimental selection criteria using CMUL 
(CM) or Erat (ER) observables. 

To check our centrality selection method we used events generated with the 
IQMD model (see Section 4), for which the impact parameter is known. Fig. 6 
shows the distribution of the impact parameters for different centralities se- 
lected in the way outlined above, for IQMD HM, Au+Au at 150-A MeV. As 
one could see, the argument for our centrality selection based on the CMUL- 
Erat correlation is confirmed by the model predictions. 



2.3 Reaction plane and flow angle determination 

Besides the collision geometry, the reaction plane and the sidewards flow an- 
gle determination are crucial for the quality of the information extracted from 
azimuthal distributions. The transverse momentum analysis method has been 
used for the reaction plane reconstruction [46]. In order to avoid the autocor- 
relations, the particle of interest was not included in the reaction plane deter- 
mination. In such a situation the momentum conservation is violated and a 
recoil correction can be done. One has to mention that such a correction brings 
back the autocorrelations and, taking into account that we will concentrate 
on heavy systems, characterized by high particle multiplicity, where the recoil 
correction is negligible, we rather preferred not to use it. 

The correction for the reaction plane fluctuations due to the finite multiplicity 
and detector biases has been estimated using the method proposed by Olli- 
trault [47]. This is based as the previous one [46] on a random subdivision 
of each event in two and calculating the difference between the azimuth an- 
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full - Ollitrault 
dashed - <cos nA<|) 12 > 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

<cos Acf> / j> 

Fig. 7. The correction factors for oi and ai (see Eq. 4 in the text) due to the reaction 
plane resolution. 

gles of the reaction planes extracted from these subevents: A0^ = $1 — $2- 
The dependence of this correction on the cosine of is represented in 

Fig. 7. In the same figure is represented the correction using (cosnA$ 12 ), 
where A$ 12 =A0#/2 is a measure of the reaction plane resolution [46]. Ta- 
ble 1 contains the (cos A<pf>) values for the energies and centralities studied in 
the present paper. 

Table 1 

Experimental values of (cosA^) for the measured energies and centrality bins. 



Beam Energy 


90 


120 


150 


250 


400 


CM2 


0.151 


0.390 


0.454 


0.756 


0.820 


CM3 


0.252 


0.520 


0.644 


0.827 


0.864 


ER4 


0.283 


0.527 


0.679 


0.797 


0.817 


ER5 


0.164 


0.326 


0.471 


0.563 


0.600 



An example of A0/j distribution for ER4 centrality is presented in Fig. 8a 
for E=250-A MeV. Such experimental representations are used for determin- 
ing the dimensionless parameter x using the following expression (Eq. 12 in 
ref. [47]): 



dN e"* /2 (2 



-(1 + X 2 /2) + *) + L (z)] + X 2 /2[Ii(z) + Mz)])(l) 

7T J 



dAfa 2 

Where z=xfcosA</># and L , Li are the modified Struve functions. The line in 
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Fig. 8a is the result of a fit using the above equation. The correction factors 
are then given by Eq. 8 in ref. [47]: 




Fig. 8. a) Reaction plane resolution for E=250-A MeV, ER4; b) (cos A</>^) depen- 
dence on the scaled multiplicity (cm5 is the lower limit for the CM5 bin; cm5=56 
for the present example). 

A representation of the (cos A0^) as a function of scaled CMUL for the energy 
of 250-A MeV is presented in Fig. 8b. Table 1 shows clearly the influence 
of the incident energy and of the impact parameter on the precision of the 
reaction plane determination. For a given incident energy the precision reaches 
a maximum for an intermediate value of the impact parameter (multiplicity) 
as it can be seen also in Fig. 8b. This figure shows also the behaviour of 
(cos A(J)r) as a function of CMUL using the information delivered by the full 
device or only by some of its subdetectors. In the subsequent analysis, the full 
detector has been used for the reaction plane reconstruction. 

Already in the first paper in which the experimental confirmation of the out- 
of-plane enhancement was reported [2] it was stressed the importance of per- 
forming the azimuthal distribution analysis in a reference system which has 
the polar axis along the sidewards flow direction [4]. There are many recipes 
to determine the flow angle. The methods based on the sphericity analysis 
[48] or on fitting the momentum distribution of the particles detected in the 
participant region by a three dimensional anisotropic Gaussian distribution 
[49] may be affected by the contributions coming from the spectator compo- 
nents. Information on flow angle can be also obtained using sidewards flow 
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analysis [52]. In the present analysis we use a different procedure for theta 
flow determination [23,50,9]. The flow angle Qfi ow is defined as being the an- 
gle with which all the events in a given centrality class have been rotated in 
order to maximize the squeeze-out signal or by carefully analyzing the p x - y 
distributions in the region of E tran value. 

Au+Au E=250 AMeV 



1.8 r 
1.6- 
1.4 
1.2- 
1 1 
0.8 
0.6 



CM3 




A 2=3 ® flow =21.8 
Z=2 ® flow =21.2 
• Z=1 @ flow =17.4 



5 10 15 20 25 30 35 40 



flow 



(deg.) 



Fig. 9. Squeeze-out ratio Rn in a reference system rotated in the reaction plane by 
"© flow" angle relative to the beam axis, integrated over all momenta. 

Fig. 9 shows the dependence of R e ^ p : 



f (0 = 90°) + ^(0 = 270°) 



rjexp 



f (0 = 0°) + ^(0 = 180° 



(3) 



estimated in a coordinate system which has two axes (x — z) in the reaction 
plane and the polar axis z at an angle u Ofi ow " relative to the beam axis, as 
function of the rotation angle "0 fi ow " for 250-A MeV incident energy and CM 3 
centrality bin. As one could see, the maximum value of Rn and the width of 
its distribution as a function of rotation angle "0/z O uj" shows a dependence on 
the mass of the analyzed particle. The mass dependence of the distribution 
width is a direct consequence of the sidewards flow profile. Part of this coming 
from thermal fluctuations, it is natural that with increasing the mass of the 
analyzed particle the distribution becomes narrower [52]. If Qfi ow extracted 
from a standard sidewards flow analysis is affected by the flow profile [52], the 
mean value of distribution as a function of u Ofi ow " does not depend on 
it. Nevertheless, the data still show a mass dependence of Qfi ow . It could be 
attributed to the dynamics of the collision. The light particles quite probable 
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are originating from the initial stage of the collision, when the sidewards flow 
is not yet developed and transparency effects are more important, leading to 
lower values of the sidewards flow angle. On the contrary, heavier fragments 
could originate from later stages of the collision, when the sidewards flow has 
developed. Part of this pattern could be also due to preequilibrium and se- 
quential emission processes. Already for Z=2 and Z=3 fragments the difference 
is negligible. Consequently, the Qfi ow values used in the present analysis cor- 
respond to these reaction products. In order to see the difference between the 
azimuthal distributions in the rotated reference system along the sidewards 
flow direction and in the non-rotated one, along the beam axis, we represent in 
Fig. 10 the corresponding two azimuthal distributions for the energy of 120-A 
MeV. We have selected here particles of mass A=4, the centrality bin ER4 and 
normalized transverse momentum per nucleon to the projectile momentum per 
nucleon in the cm. system p[ ^ = (p t / A)/(p c ™/ A P ) >0.8. 



The effect of the rotated reference frame on the Rn values makes unnecessary 
any other argument for the need of carrying out these types of studies in this 
reference system. 

Au+Au E=120 AMeV ER4 A=4 



=3 




o 



• flow axis ( o reflected ) 

* beam axis 



pM >0.8 



90 180 270 360 

4 (deg.) 



Fig. 10. Azimuthal distributions for A=4 particles in the beam and in the rotated 
reference frame, 80° < G cm <100°. 



2.4 The influence of the detector acceptance 



The FOPI device having an azimuthal symmetry [36,38,39], the border re- 
gions between different subdetectors are symmetric relative to the beam axis. 
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In other words, the transition between two consecutive subdetectors, which 
causes small shadows (see Fig. reffig-1), takes place for a given rapidity at a 
given transverse momentum, and remains constant as a function of azimuth 
(4>). Once the reference frame is rotated, the symmetry is broken and the <fi 
dependence of these regions have to be carefully treated. Fig. 11 shows as 
an example how the 0^=34° region (at midrapidity in the rotated system, 
pl=0) which corresponds to the borderline between CDC and forward Plastic 
Wall is seen in a — <fi representation in reference frames rotated relative to 
the beam axis by different Qfi ow values. 

& 2 
Ct 1.8 

1.6 

1.4 

1.2 
1 

0.8 
0.6 
0.4 
0.2 



"0 90 180 270 360 

<t> (deg.) 

Fig. 11. The geometrical thresholds as function of the azimuthal angle in the rotated 
reference frame for different flow angle values. 

Due to this reason, for pf^ >0.8, where the phase space is covered mainly 
by the CDC, the analyzed azimuthal range was [0°,90°] and [270°,360°]. The 
complete distributions have been obtained by reflecting these ranges relative 
to 90° and 270°, respectively. 



3 Experimental Results 

In the present paper the azimuthal distributions are analyzed in a reference 
frame rotated in the reaction plane by Qfi ow - The midrapidity region is selected 
by the polar angular range in the rotated reference frame 0{ m =8O°-lOO°. For 
the present study, concentrated on the transition energy such a selection is 
more adequate than one based on p z or rapidity. 
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3.1 Flow angle 



The incident energy dependence of fi ow determined using the recipe described 
in the previous chapter is shown in Fig. 12 for different centralities. The source 
of the slight differences relative to the ones obtained using transverse momen- 
tum analysis [52] has been discussed in the previous chapter. An increase of 
Qfiow as function of energy and centrality up to 400- A MeV and ER4, respec- 
tively, can be followed in Fig. 12. The lines, drawn in order to guide the eyes, 
represent the result of polynomial fits to the experimental points. 




E inc (MeV/u) 

Fig. 12. The excitation function of the flow angle for the four centrality bins. 

As far as 90- A MeV is in the range of E tran values (see later), the determination 
of Qfi ow using the present procedure is not very accurate. While this does not 
influence the results relative to E tran studies, it could bring a big uncertainty 
in the determination (by extrapolation) of the balance energy, Ej, a i. Although 
this subject will be properly treated in a separate paper, it is worth to mention 
that the present analysis gives larger E ba i values than those reported in ref. 
[58], but in agreement with our previous results [12]. 

3.2 General trends of the azimuthal distributions 

For a multidimensional study of the incident energy at which a transition 
from an in-plane enhancement to the out-of-plane preferential emission takes 
place, the excitation functions of azimuthal distributions for different frag- 
ment species, transverse momenta and centrality have been analyzed. We will 
present some of these trends in what follows. 



15 





40 






-2. 


35 






f 






30 






25 




20 




15 




10 




5 








Au+Au E=150 AMeV CM3 



Au+Au E=150 AMeV ER4 



A=2 



1 .4<p,(0) <1 .6 
1 .2<p,(0) <1 .4 
1 .0<p t (0) <1 .2 
0.8<p,(°) <1 .0 

li 





reflected * 



270 360 

<|> (deg.) 

Fig. 13. The azimuthal distributions for 
different ranges of the scaled momentum 
pf^ for 150-A MeV incident energy and 
A=2 reaction products. 
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Fig. 14. The azimuthal distributions as 
function of the mass of the detected par- 
ticle for 150-A MeV and ER4 centrality. 



For a given incident energy (as an example 150-A MeV was taken), the de- 
pendence of the azimuthal distributions as a function of the scaled transverse 
momentum pf 1 ' is presented in Fig. 13 for CM3 centrality and A=2 reaction 
products. The enhancement of the anisotropy with increasing transverse mo- 
mentum seen in the previous experiments [21,22,11] is confirmed. Taking the 
advantage of a complete phase space coverage of the present configuration of 
our experimental device, we chose to do the present studies for high transverse 
momentum p[°\ a region where the theoretical predictions [26,30] suggested 
that the squeeze-out signal is sensitive to the nuclear matter Equation of State 
(EoS). For a given incident energy and range of pf^ the dependence of the az- 
imuthal distribution as function of mass of the analyzed particle is shown 
in Fig. 14. The well known enhancement of the azimuthal anisotropy with in- 
creasing the mass of the analyzed reaction product is evidenced [4,26,22,23,11]. 

Fig. 15 shows once more the importance of doing these types of studies in the 
appropriate coordinate frame [2]. In the laboratory system the dependence 
of the squeeze-out signal as function of the centrality is biased as far as one 
looks to the projection of these distributions in a plane which is not normal to 
the flow direction. As one can see from Fig. 12 the flow angle increases with 
the centrality up to ER4 and this influences the information on the squeeze- 
out signal as a function of centrality if this information is extracted from the 
reference frame with the polar axis along the beam direction. This caution has 
to be taken when one studies the excitation function of the squeeze-out signal. 
A sample of the incident energy dependence of the azimuthal distributions for 
A=4 reaction products and for a given range in transverse momentum, 0.8< 
pf^ <1.8, is presented in Fig. 16. All the azimuthal distributions presented are 
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Fig. 15. Centrality dependence of the 
azimuthal distributions for 150-A MeV, 



A=4 fragments and 0.8< 
transverse momentum. 



rf» <1.8 
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Fig. 16. Energy dependence of the az- 
imuthal distributions for ER4 central- 



ity, A=4 fragments and 0.8< p\ 
transverse momentum. 



(0) 



<1.8 



the primordial ones, namely not corrected for the reaction plane resolution. 

As it was already mentioned, the azimuthal distributions have been analyzed 
in the rotated reference frame. They have been fitted with a second order 
Fourier expansion: 

fdN\ exp 

—— = Orf p ■ (1 + a{ xp cos + a^ p cos 20) (4) 
\d(p ) 



The smooth curves in Fig. 13-16 have been obtained using this fit function. As 
there are some systematic dips around 90o and 270o, caused by the detector 
acceptance in the rotated reference system and by the influence of the magnetic 
field on the reaction plane reconstruction, we excluded them from the fits, a® 
gives the average value of the distribution, a\ corresponds to in-plane flow 
and ci2 describes the squeeze-out pattern. The squeeze-out ratio R e ^ p can be 
written: 

-r exp 

Tjexp _ ~ a 2 /r\ 
~ i , exp \° 

1 + a 9 



The a^ v values have been corrected for the fluctuations of the estimated re- 
action plane using the method of ref. [47] explained in the previous chapter: 

a -- = a -7(cos2A0> (6) 
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The corrected quantity Rn, which will be used in the following, is given by: 

i n corr 

Rn = 1 a * (7) 



£ 02 

I 
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S 

to 



Au+Au Z=1+2 



C 



FOPI CM3 



• r.p. corr. 

o no r.p. corr. 

* Plastic Ball 



in-plane 

out-of-plane 




200 400 600 800 1000 

Beam energy (MeV/u) 

Fig. 17. The excitation function for the 02 coefficient for the CM3 centrality (see 
text for details). 

Before proceeding to the detailed analysis of the transition energy, we present 
in Fig. 17 the excitation function for the a2 coefficient extracted using Z=l and 
2 particles and integrated over the transverse momentum. To compare with 
the Plastic Ball results [4] which were not corrected for the reaction plane 
fluctuations, our a 2 values are shown both with (full dots) and without (open 
circles) this correction. Both sets of data are for roughly similar centrality bin 
(CM3 for the present analysis and MUL3 for the Plastic Ball data). Only for 
this case, we used a p z window to select the midrapidity, as it was done in 
ref. [4], namely \{p{ / A) / {p^ / A P )\ <0.1, where p{,/A is the particle longitudi- 
nal cm. momentum per nucleon in the rotated reference frame and p c ™/Ap is 
the projectile momentum per nucleon in the cm. system. A good agreement 
between the two data sets can be observed. This agreement is meaningful as 
the resolutions of the reaction plane for the two experiments are comparable, 
at least at 400-A MeV [4]. 

The experimental a£ v values in the rotated reference frame at the incident 
energies, centralities, transverse momenta and for reaction species studied in 
this paper are given in the Appendix. 
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3.3 Multidimensional analysis of the transition energy 

Most of the previous studies of the azimuthal distributions have been con- 
centrated on the behaviour of the observed anisotropies as a function of cen- 
trality, rapidity, transverse momentum and the mass of the reaction products. 
A detailed study of the incident energy at which a transition from in-plane 
to out-of-plane emission takes place could give deeper insight on the relative 
contribution of the attractive and repulsive forces, the lifetime of the emitting 
source, its rotational energy and expansion dynamics. Earlier evidences of such 
a transition energy have been reported for Au+Au combination [6,9,11,12] and 
Zn+Ni [5]. These studies have been done in the reference frame in which the 
polar axis was along the beam direction. At 100-A MeV Au+Au collisions 
a transition from squeeze-out at high centrality to a rotational-like pattern 
towards peripheral collisions was observed [8]. 

2 3 
DC 

2.5 
2 
1.5 
1 



v " 100 200 300 400 

E inc (MeV/u) 

Fig. 18. The excitation function for the squeeze-out ratio for CM3 centrality and 
transverse momentum window 0.8< <1.0. 

A multidimensional analysis of this transition energy in terms of its depen- 
dence as function of centrality, transverse momentum and mass of the reaction 
products for Au+Au is performed in this section. The procedure used to ex- 
tract the transition energy (E tran ) can be followed in Fig. 18. For a given 
centrality range, transverse momentum and mass of the reaction product, the 
Rn value is represented as a function of incident energy (E inc ). Squeeze-out 
signal corresponds to Rn > 1 while in-plane enhancement of the azimuthal 
distribution is characterized by Rn < 1. The continuous lines correspond to 
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the result of fits of the experimental points for different masses using a poly- 
nomial function. The intersections of these lines with the one corresponding 
to Rn = 1 value are defined as E tran . Fig. 18 shows a clear mass dependence 
of the transition energy, lower values corresponding to heavier particles for a 
given centrality and a given range in the scaled transverse momentum p[ . 
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Fig. 19. The transition energy as a function of the scaled momentum for different 
particle types for the CM2 and CM3 centrality bins, with (left column) or without 
(right column) reaction plane correction. 

Following the above recipe for different regions of pf\ a two dimensional 
dependence of the E tran values as function of mass of the analyzed reaction 
products and transverse momentum can be obtained. The results are presented 
in Fig. 19 for CM2 and CM3 centrality bins, with and without doing the 
correction for the reaction plane resolution. As expected, in case of using the 
correction, the values of E tran are systematically slightly higher, but the trends 
are identical for both cases. A continuous decrease of the E tran as function of 
Pt is evidenced for all analyzed particles and the difference in E tran values for 
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different particles is decreasing towards larger values of p t ' for both centrality 
bins. One should mention here that for the cases when < 1 values have not 
been reached E tran was determined by extrapolating the fits of as function 
of incident energy (E inc ). 

For a rotating emitting source one would expect a larger in-plane alignment for 
heavier fragments [7]. This effect alone can not explain the mass dependence 
of E tran . Therefore a dynamical effect has to be considered besides the pure 
geometrical one of shadowing. Different particles, originating from different 
regions of the fireball [34,53] would feel the shadowing in a different way. 
At large p[ ^ the contribution comes from larger expansion velocities, earlier 
expansion phase of the fireball [34,38] and consequently higher shadowing. At 
lower values of pf \ the expansion zones are less localized, especially for the 
light fragments due to larger contribution of the thermal velocities relative to 
the collective ones, the shadowing being less effective on the in-plane yields. 
Contaminations coming from sequential evaporation processes smear-out such 
effects but cannot explain the observed trends. 



4 IQMD calculations 

The IQMD (Isospin Quantum Molecular Dynamics) [44] is a QMD model [25] 
which takes into account the isospin degree of freedom for nucleon-nucleon 
cross section and Coulomb interaction. The main argument for using this ver- 
sion in the present context is based on the fact that IQMD has been success- 
fully used for analyzing flow phenomena in heavy ion collisions [26,27,30,56]. 
In this work two different parametrizations of the EoS are used, a hard EoS 
(compressibility K— 380 MeV) with momentum dependence of the nucleon 
interaction (MDI) - HM and a soft EoS (K= 200 MeV), with MDI - SM. 

The events produced by the model have been filtered by the experimental 
filter and analyzed in a similar way as the experimental data. Table 2 is the 
analog of Table 1 but this time representing the resolution of the reaction plane 
determination for the events generated by this model (HM parametrization). 
The worse resolution at lower energies for calculated events is the natural 
consequence of lower yield for intermediate mass fragments predicted by the 
model. This was the reason why for the calculations we used the true reaction 
plane. 

The model was used to estimate the accuracy of the reaction plane corrections 
[46,47] and it has been found that the correction used in the present study 
works fine in what concerns the sidewards flow, but it is overestimated for 
the squeeze-out signal, especially for the cases (cos A0#) <0.4. We also found 
that applying the correction in the rotated frame slightly overestimates the 
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Table 2 

The values of (cos A^r) obtained analyzing the IQMD (HM) model events. 



Beam Energy 


90 


120 


150 


250 


400 


CM2 


0.005 


0.087 


0.378 


0.771 


0.859 


CM3 


0.036 


0.309 


0.576 


0.837 


0.887 


ER4 


0.081 


0.395 


0.579 


0.800 


0.863 


ER5 


0.074 


0.161 


0.266 


0.498 


0.543 



real squeeze-out signal. 
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Fig. 20. The excitation function for the squeeze-out signal for CM3 centrality bins, 
data and IQMD model for different particle selections. 

As mentioned already, it is well known that the IQMD model underpredicts 
the composite particles yield [57,35] and this drawback is especially important 
for our energy range. To establish the most meaningful way of comparison, in 
Fig. 20 we compare the excitation function of the squeeze-out ratio (Rn) for 
data and the model for three ways of particle selection: a) protons only, b) 
all particles up to alphas and c) a coalescence- type, summing and weighting 
the azimuthal distribution by each particle's atomic number A (including up 
to alphas here as well). It was shown within the IQMD model [30] that the 
sensitivity of the squeeze-out signal to EoS is present only for high transverse 
momenta. With our selection, 0.8 < < 1.8, at the higher energies we do see 
a dependence of the squeeze-out on the EoS. Although the model predicts an 
increase of the squeeze-out signal with particle mass [26] , because the protons 
are dominating in IQMD the model values for the three cases show only small 
variations. This is clearly not the case with the data, where the inclusion of 
the composites, without or with weighting, is increasing the squeeze-out sig- 
nal significantly. For protons the model predicts larger squeeze-out than the 
experimental values which compensates the lower yield of composite particles 
produced by a coalescence mechanism, so that integrating over A<4 particles 
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the agreement with the experiment is good. However, Fig. 20 shows that the 
nucleonic squeeze-out type of flow predicted by the model and the coalescence 
mechanism used to produce fragments do not explain the experimental ob- 
servations. This conclusion holds for the ER4 centrality as well. In Fig. 21 
we present for this centrality the comparison for A<4 particles weighted by 
their mass. Apart of different magnitudes, the behavior is similar to the one 
presented in Fig. 20c. All these facts show the importance of a comparison be- 
tween the experiment and model either for each particle or for coalescence- type 
(as in Fig. 20c and Fig. 21) and not simply summing-up all particles, when 
different effects could compensate each other and the agreement is artificially 
good. 

Au+Au ER4 0.8 < p,(0) < 1.8 

DC 

2 
1.5 
1 



w " 100 200 300 400 

E inc (MeV/u) 

Fig. 21. The excitation function for the squeeze-out signal for ER4 centrality bin, 
data and IQMD model for A<4 particles weighted by their mass. 

Although at low energies the compressibility should indeed not play an impor- 
tant role, there is a systematic difference between the two parametrizations, 
the HM version predicting a higher in-plane alignment. Similar trends have 
been observed at 100-A MeV for Au+Au by Aladin-Miniball Collaboration 
and comparisons with model predictions using BUU calculations show bet- 
ter agreement with a stiff EoS, momentum dependent nuclear interaction and 
reduced in-medium cross section (<Jnn = 0.8 x o"/ ree ) [8]. The in-plane en- 
hancement in Ar+V system at energies below 100-A MeV was studied using 
BUU calculations with a nn = 0.5 x <T/ ree and a stiff EoS [7]. More recently the 
transition region was investigated for Ca+Ca system using an isospin depen- 
dent BUU model and the dependences of E tran on EoS, MDI and ctjvjv were 
found to be significant [14]. 

Although the number of generated events was about 2000 per 1 fm interval 
in the impact parameter, the statistics for particles heavier than A=l was 
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Fig. 22. Experimental transition energy values as function of centrality for: a) 
Z=l,2,3 particles and b) A=l,2,3,4 particles. Full symbols denote the experimental 
values, open symbols are the IQMD model predictions: crosses - HM and stars - SM 
parametrization for the EoS. 

not sufficient to perform with reasonable accuracy a similar analysis of E tra n 
as the one done for the experimental data. Consequently, we have chosen to 
do the comparison in a broad pt interval and as function of centrality. The 
centrality dependence of the E tran values is presented in Fig. 22, for 0.8 < 
p[ ^ < 1.8 window and particles selected according to: a) charge (Z=l-3) and 
b) mass (A=l-4). For Z=l and A=l, respectively, the experimental data are 
compared with the IQMD model with the two versions for compressibility. The 
shown centrality bins (ER5, ER4, CM3 and CM2) correspond to the impact 
parameter intervals 0-2, 2-4, 4-6 and 6-8 fm, respectively, as it was explained in 
the second chapter of this paper. Systematically, the SM version predicts lower 
transition energies relative to HM version. The model predicts an evolution of 
the E tran as function of centrality for Z=l and A=l particles similar to the 
experiment. With decreasing impact parameter the amount of energy which 
can be transferred in a rotational-like motion becomes smaller. At the same 
time the volume of the spectator matter causing the shadowing is reduced, 
while the number of baryons involved in the fireball region, A part , increases, 
contributing to larger expansion [51,54,55,37]. All these effects play a role in 
the observed experimental azimuthal distributions. However, the centrality 
dependence of the E tran suggests that the reduction of the rotational motion 
of the fireball is the main source of the observed effect. One could observe in 
Fig. 22 a less pronounced dependence on centrality for Z=l particles relative to 
the heavier ones. A possible explanation would be that at higher centralities, 
the three sources of emitting particles, the projectile spectator, the target 
spectator and the participant zone, are less separated, the mixed contribution 
being enhanced for lighter reaction products. The contribution from sequential 
processes brings an additional smearing of anisotropies, increasing artificially 
the E tran value for lighter particles. 
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5 Conclusions 



Using a new generation of data collected with the complete FOPI experimen- 
tal device we presented a systematic study of the excitation function of the 
azimuthal distributions around midrapidity. For Au+Au system for which the 
most complete excitation function was measured, starting with 90-A MeV, a 
clear transition from in-plane (rotational-like) to out-of-plane (squeeze-out) 
was evidenced to take place going from lower to higher energies. The inci- 
dent energy at which this transition takes place, called E tran , was studied 
as function of mass of reaction products, their scaled transverse momentum 
(Pt°^) an d centrality. Special attention was payed to the sidewards flow angle 
inflow) determination used to define the rotated reference frame in which all 
these studies have been done. 

In general one could say that a systematic decrease of the E tran values with 
increasing the mass of the analyzed fragment, the transverse momentum and 
centrality of the collision was observed. E tran values being the incident en- 
ergies where the effect of the rotation of a compound system created by the 
mean field, expansion of a hot and compressed participant zone built-up by 
many consecutive nucleon-nucleon interactions and the shadowing effect of the 
colder spectator regions compensate each other, it seems to be an experimen- 
tal observable which could distinguish between different theoretical models 
and different approximations used by them. Studying an expanding-rotating 
hot and compressed fireball in the presence of the shadowing objects can give 
additional information on the expansion dynamics. The observed decrease of 
the Etran as a function of the mass of the emitted fragments can be the result 
of a complex dynamics, indicating that heavier particles for which the thermal 
contribution is less important are better localized in the expansion process us- 
ing selection than the lighter ones. When integrating over a large range 
in momenta, our E tran values are in a good agreement with previous studies 
which found the transition to take place around 100-A MeV. 

Detailed comparisons with the IQMD model predictions evidenced the ne- 
cessity to avoid any kind of fragment mixture in order to draw significant 
conclusions. Squeeze-out ratio for protons and coalescence-type revealed that 
nucleonic out-of-plane flow and the coalescence mechanism used by the model 
to create complex fragments can not explain the experimental observations. 
The low statistics of IQMD generated events did not allow a similar analysis 
of E tran as for the experiment. Integrating on pf\ the systematic decrease 
of Etran values with increasing the collision centrality is reproduced by the 
model. Although systematically the hard EoS gives larger E tra n values, closer 
to the experimental values, the small difference relative to the soft EoS is not 
sufficient to draw definite conclusions on the EoS. 
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Appendix: Experimental a 2 values 

In this Appendix we present the values of the a^ p coefficient for different 
particle types and ranges in normalized transverse momenta. Remind that 
these values are the ones obtained from the fit of the experimental data with 
expression 4, namely without including the correction for the reaction plane 
resolution. In parentheses we quote the absolute errors on the last digit of the 
respective number. 

Table 3 

The a^ p values in the rotated reference frame for particles with A=l, 2 and 3, for 
CM2 centrality bin. 



(0) 

Pt 


A 


90 


120 


150 


250 


400 


0.8-1.0 


1 


0.034(9) 


0.044(9) 


0.020(9) 


-0.030(6) 


-0.110(6) 




2 


0.07(1) 


0.06(1) 


-0.02(1) 


-0.230(7) 


-0.400(8) 




3 


0.08(1) 


0.01(1) 


-0.12(1) 


-0.48(1) 


-0.71(1) 


1.0-1.2 


1 


0.047(8) 


0.030(9) 


-0.004(9) 


-0.110(6) 


-0.210(7) 




2 


0.06(1) 


0.00(1) 


-0.12(1) 


-0.400(9) 


-0.59(1) 




3 


0.05(2) 


-0.09(2) 


-0.31(2) 


-0.71(2) 


-0.87(2) 


1.2-1.4 


1 


0.044(9) 


0.013(9) 


-0.060(9) 


-0.180(7) 


-0.310(8) 




2 


0.04(1) 


-0.07(2) 


-0.24(2) 


-0.55(1) 


-0.71(2) 




3 


0.03(2) 


-0.27(4) 


-0.55(5) 


-0.88(2) 


-0.81(6) 


1.4-1.6 


1 


0.03(1) 


-0.050(9) 


-0.12(1) 


-0.310(8) 


-0.44(1) 




2 


0.01(2) 


-0.24(3) 


-0.40(3) 


-0.72(2) 


-0.89(3) 
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Table 4 

Same as Table 3, but for CM3 centrality bin. 



(0) 

Pt 


A 


90 


120 


150 


250 


400 


0.8-1.0 


1 


0.018(9) 


0.022(9) 


0.007(9) 


-0.030(6) 


-0.060(6) 




2 


0.04(1) 


0.04(1) 


-0.04(1) 


-0.200(6) 


-0.310(7) 




3 


0.06(1) 


-0.04(1) 


-0.18(1) 


-0.45(1) 


-0.62(1) 


1.0-1.2 


1 


0.021(8) 


0.025(9) 


-0.001(9) 


-0.070(6) 


-0.140(6) 




2 


0.05(1) 


-0.01(1) 


-0.12(1) 


-0.340(7) 


-0.470(9) 




3 


-0.00(2) 


-0.15(2) 


-0.35(2) 


-0.63(2) 


-0.77(2) 


1.2-1.4 


1 


0.023(8) 


-0.010(9) 


-0.060(9) 


-0.140(7) 


-0.230(8) 




2 


0.01(1) 


-0.11(2) 


-0.26(2) 


-0.48(1) 


-0.60(1) 




3 


-0.01(2) 


-0.29(3) 


-0.51(3) 


-0.78(3) 


-0.89(3) 


1.4-1.6 


1 


0.014(9) 


-0.040(9) 


-0.14(1) 


-0.230(8) 


-0.35(1) 




2 


-0.03(2) 


-0.20(2) 


-0.39(2) 


-0.62(2) 


-0.71(2) 




3 


-0.10(5) 


-0.49(6) 


-0.60(6) 


-0.8(2) 





Table 5 

Same as Table 3, but for ER4 centrality bin. 



(0) 

Pt 


A 


90 


120 


150 


250 


400 


0.8-1.0 


1 


0.021(9) 


0.034(9) 


0.012(9) 


-0.026(6) 


-0.068(6) 




2 


0.04(1) 


0.04(1) 


-0.05(1) 


-0.19(1) 


-0.24(1) 




3 


0.04(1) 


-0.04(1) 


-0.12(1) 


-0.32(1) 


-0.48(1) 


1.0-1.2 


1 


0.019(8) 


0.037(9) 


0.008(9) 


-0.068(6) 


-0.12(1) 




2 


0.04(1) 


-0.02(1) 


-0.09(1) 


-0.24(1) 


-0.37(1) 




3 


-0.02(2) 


-0.10(2) 


-0.30(2) 


-0.49(1) 


-0.61(2) 


1.2-1.4 


1 


0.014(8) 


0.018(9) 


-0.052(9) 


-0.13(1) 


-0.21(1) 




2 


0.02(2) 


-0.06(2) 


-0.17(2) 


-0.42(1) 


-0.49(1) 




3 


-0.05(3) 


-0.22(3) 


-0.35(3) 


-0.66(2) 


-0.75(2) 


1.4-1.6 


1 


0.029(9) 


-0.074(9) 


-0.09(1) 


-0.20(1) 


-0.26(1) 




2 


-0.04(2) 


-0.20(2) 


-0.25(2) 


-0.52(2) 


-0.58(2) 
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Table 6 

Same as Table 3, but for ER5 centrality bin. 



(0) 

Pt 


A 


90 


120 


150 


250 


400 


0.8-1.0 


1 


0.03(1) 


0.02(2) 


0.00(2) 


-0.01(2) 


-0.03(1) 




2 


0.03(1) 


0.03(1) 


-0.01(2) 


-0.09(2) 


-0.15(1) 




3 


0.02(1) 


0.00(2) 


-0.05(2) 


-0.11(1) 


-0.23(1) 


1.0-1.2 


1 


0.00(1) 


0.01(1) 


-0.01(2) 


-0.04(2) 


-0.05(2) 




2 


0.01(1) 


0.00(1) 


-0.04(2) 


-0.12(2) 


-0.21(2) 




3 


0.00(2) 


-0.04(2) 


-0.13(2) 


-0.25(2) 


-0.34(2) 


1.2-1.4 


1 


0.01(1) 


-0.03(2) 


-0.05(2) 


-0.08(2) 


-0.11(2) 




2 


0.01(2) 


-0.03(2) 


-0.11(2) 


-0.19(2) 


-0.28(2) 




3 


0.01(2) 


-0.08(2) 


-0.14(2) 


-0.36(2) 


-0.39(4) 



Table 7 

The a^ v values in the rotated reference frame for particles with Z=l, 2 and 3, 
integrated over the momentum interval 0.8< p{°^ <1.8. 



Centrality 


Z 


90 


120 


150 


250 


400 


CM2 


1 


0.040(4) 


-0.022(4) 


-0.084(4) 


-0.240(3) 


-0.330(3) 




2 


0.02(1) 


-0.23(2) 


-0.43(2) 


-0.82(2) 


-0.93(2) 




3 


0.01(1) 


-0.37(4) 


-0.64(5) 


-0.84(4) 




CM3 


1 


0.022(4) 


-0.034(3) 


-0.120(3) 


-0.241(2) 


-0.298(3) 




2 


-0.04(1) 


-0.33(1) 


-0.56(1) 


-0.80(1) 


-0.89(1) 




3 


-0.06(2) 


-0.43(2) 


-0.72(3) 


-0.97(3) 




ER4 


1 


0.017(6) 


-0.029(5) 


-0.10(4) 


-0.190(3) 


-0.240(3) 




2 


-0.05(1) 


-0.29(1) 


-0.49(1) 


-0.68(1) 


-0.78(1) 




3 


-0.09(3) 


-0.40(2) 


-0.61(2) 


-0.84(2) 


-0.88(4) 


ER5 


1 


0.008(6) 


-0.010(9) 


-0.05(1) 


-0.10(1) 


-0.14(1) 




2 


-0.02(1) 


-0.15(2) 


-0.29(2) 


-0.43(2) 


-0.54(2) 




3 


-0.06(2) 


-0.18(5) 


-0.36(4) 


-0.48(4) 


-0.63(5) 
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Table 8 

The a^ p values in the rotated reference frame for particles with A=l,2,3 and 4, 

integrated over the momentum interval 0.8< p{°^ <1.8. 

Centrality A 90 120 150 250 400 

CM2 1 0.034(6) -0.001(4) -0.050(4) -0.170(6) -0.260(6) 

2 0.052(6) -0.010(6) -0.110(6) -0.361(4) -0.50(1) 

3 0.06(1) -0.05(1) -0.19(1) -0.55(1) -0.74(2) 

4 0.02(1) -0.23(2) -0.43(2) -0.82(2) -0.93(3) 
CM3 1 0.017(6) -0.013(4) -0.070(8) -0.160(4) -0.221(4) 

2 0.025(6) -0.05(1) -0.15(1) -0.34(1) -0.43(1) 

3 0.03(1) -0.10(1) -0.27(2) -0.52(1) -0.66(1) 

4 -0.04(1) -0.33(2) -0.56(2) -0.80(2) -0.89(2) 
ER4 1 0.018(8) -0.007(9) -0.061(8) -0.110(6) -0.160(6) 

2 0.023(8) -0.03(1) -0.11(1) -0.27(1) -0.33(1) 

3 0.00(1) -0.06(1) -0.20(2) -0.39(1) -0.52(1) 
3 -0.05(1) -0.29(2) -0.49(2) -0.68(2) -0.78(2) 

ER5 1 0.006(6) -0.005(9) -0.027(9) -0.044(8) -0.079(8) 

2 0.011(8) -0.03(1) -0.05(1) -0.13(1) -0.20(1) 

3 0.01(1) -0.04(2) -0.09(2) -0.17(1) -0.29(2) 

4 -0.02(1) -0.15(1) -0.29(1) -0.43(1) -0.54(2) 
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